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Theories for new particle and early-Universe physics abound with pseudo-Nambu-Goldstone fields 
that arise when global symmetries are spontaneously broken. The coupling of these fields to 
the Chern-Simons term of electromagnetism may give rise to cosmological birefringence (CB), a 
frequency-independent rotation of the linear polarization of photons as they propagate over cosmo- 
logical distances. Inhomogeneities in the CB-inducing field may yield a rotation angle that varies 
across the sky. Here we note that such a spatially-varying birefringence may be correlated with 
the cosmic microwave background (CMB) temperature. We describe quintessence scenarios where 
this cross-correlation exists and other scenarios where the scalar field is simply a massless specta- 
tor field, in which case the cross-correlation does not exist. We discuss how the cross-correlation 
between CB-rotation angle and CMB temperature may be measured with CMB polarization. This 
measurement may improve the sensitivity to the CB signal, and it can help discriminate between 
different models of CB. 
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I. INTRODUCTION 

Much attention has focused recently on cosmological 
birefringence (CB), a frequency- independent rotation of 
the linear polarization of a photon that propagates over 
cosmological distances [l| . The rotation may arise if the 
pseudoscalar of electromagnetism F'^'^F^i^ is coupled to 
a pseudo-Nambu-Goldstone field (PNGB field) that has 
variations on cosmological distances or timescales. This 
field may be identified with the quintessence fieldjl, Q 
introduced to account for cosmic acceleration 0-01 • In 
fact, the flatness required of the quintessence potential 
is naturally accommodated if quintessence is a PNGB 
field 0, 111- However, the scalar field may have nothing 
to do with quintessence — any PNGB field is expected to 
have such a coupling |9| . There may also be dark-matter 
mechanisms for CB |1C| . 

Cosmological birefringence has been sought with po- 
larized cosmological radio sources P, [ll|, [l^l , but here we 
focus on cosmic-microwave-background (CMB) probes of 
CB [l^ . If the CB-rotation angle a is uniform across the 
sky, as may result from the homogeneous evolution of 
quintessence, then there are parity- violating EB and TB 
correlations between the CMB temperature (T) and the 
curl-free (E) and curl (B) components of the CMB polar- 
ization |l3ll. Such a rotation has been sought for several 
years [ij] , and the tightest current limits on the rotation 
angle, -1.41° < a < 0.91° (95% C.L.), come from a com- 
bined analysis of the WMAPjp, BICEP [H, [13, and 
QUaD [H ^ experiments [l|T^ It is worth noting, in 
the context of the current constraints, that the uniform 
rotation angle is generally nonzero in quintessence mod- 
els for CB, while the massless-scalar- field models have no 
homogeneous time evolution and thus predict no uniform 
rotation. 

It has been pointed out that 0, more gener- 

ally, the CB-rotation might be anisotropic, giving rise to 
a(n) as a function of direction n in the sky. Refs. ^ 



|22| showed how measurement of the characteristic non- 
Gaussianities in the CMB polarization induced by a 
spatially-varying CB can be used to reconstruct a(n) 
from the CMB. The current best constraint to the root- 
variance of a, {iAa)^) ' < 4° , comes from observations 
of active galactic nuclei ■ 

In this paper we explore the possibility that CB may 
be correlated with primordial density perturbations and 
thus also with temperature fiuctuations in the CMB. 
Such correlations are to be expected, for example, if the 
CB-inducing field is a quintessence field with adiabatic 
primordial perturbations seeded during infiation. On the 
other hand, correlations between the CB angle and pri- 
mordial perturbations may be absent if, for example, the 
CB-inducing field is a massless scalar 

We first work out the predictions for spatially- varying 
a(n) for a massless-scalar- field model in which there is 
no uniform rotation. In this case, the rotation-angle 
pattern a(n) is completely uncorrelated with the CMB- 
temperature pattern T(n), and so we calculate only the 
CB-angle autocorrelation power spectrum C£" . We then 
move on to quintessence models in which the aT cross- 
correlation exists and calculate this cross-correlation 
power spectrum Cf^ . 

We derive the minimum-variance estimators for the 
aT cross-correlation and estimate the detectability of 
CB-angle fluctuations and CB-angle-temperature corre- 
lations with current and forthcoming CMB experiments. 
We find that the cross-correlation can help improve the 
sensitivity of experiments to a signal in some cases where 
the signal would otherwise be only marginally detectable. 
We show that experiments like SPIDER [2^ and Planck 
[2^ may be able to detect a cross-correlation if the CB 
signal is near its current upper bounds, while the cross- 
correlation may be detectable with a future experiment, 
like CMBPol/EPIC even if the CB power spectrum 
is several orders of magnitude smaller than the current 
upper limit. 
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This paper is organized as follows: In Section |IT] we 
introduce the two anisotropic-CB scenarios and their pa- 
rameters and calculate the corresponding cxa and aT 
power spectra. In Section IIIII we discuss how a can 
be reconstructed from a CMB temperature/polarization 
map, and then how the cross-correlations can be mea- 
sured. Here, we present the expressions for the minimum- 
variance estimators of the auto-correlation and cross- 
correlation power spectra, and expressions for their vari- 
ances. We then evaluate those variances for SPIDER, 
Planck, and CMBPol/EPIC and estimate detectability 
thresholds for these three experiments. A summary and 
concluding remarks are presented in Section IIVI Ap- 
pendix |A] details the evolution of the scalar- field pertur- 
bations, and Appendix |B] provides the full expression for 
the variance of the aT cross-correlation. 



II. SCENARIOS FOR ANISOTROPIC 
ROTATION 



We consider theories of a cosmic scalar 



coupled 



to the Chern-Simons term of electromagnetism via the 
Lagrangian 



4 2M • 



(1) 



where F, 



PP'^ /2 is the dual of the electromag- 
netic tensor, e^^pcr is the Levi-Civita tensor (totally an- 
tisymmetric), and M is a parameter with dimensions of 
mass. If is a PNGB field, then M is the vacuum ex- 
pectation value for the broken global symmetry, and /3 
is a coupling 0, @] . Such a parity-violating term in the 
Lagrangian introduces a modification of Maxwell's equa- 
tions that results in different dispersion relations for left- 
and right-circularly polarized photons. Consequently, 
linearly-polarized electromagnetic waves that propagate 
over cosmological distances undergo CB, a frequency- 
independent rotation of the plane of polarization by an 
angle a, where [l| 

a = — dr I n • 

M \dT 



£ 

M 



(2) 



where Acf) is the change in over the photon trajectory, 
and T is the conformal time. For the CMB, the polar- 
ization rotation is determined by the change in since 
recombination, when the CMB polarization pattern was 
largely established. 

Allowing for spatial fluctuations 6(f> in the cosmic scalar 
field, the anisotropy in the CB-rotation angle is then 
Aa(n) = (/3/M)(50(n), evaluated at recombination. 

Below we consider two scenarios for the scalar field. 
In the first, the scalar is massless, with no homoge- 
neous time evolution, while in the second, the scalar is 



quintessence. In both cases, CB-angle fluctuations arise 
from scalar-fleld fluctuations at the surface of last scatter 
(LSS). 



A. Massless scalar field 

In the flrst scenario, we suppose that the (j) field is 
simply a massless scalar with a potential that vanishes, 
V = 0. In this case, the value of the field is completely 
uncorrelated with primordial density perturbations^ [9|. 
If (p is eff'ectively massless during inflation there will be 
a scale-invariant power spectrum of perturbations to (j), 
Ps4,{k) ~ H] /2k^, with an amplitude flxed by the Hubble 
parameter Hj evaluated during inflation^. If we split the 
field into a smooth background component and a per- 
turbation on top of it, 4>(t,x) ~ 4'oiT) + 6(p{T,x), the 
evolution of the homogeneous component is given by the 
following equation of motion 



00 + 27^.^0 + a V - 0, 



(3) 



where T-L = a/ a, a is the scale factor, and dots denote 
derivatives with respect to conformal time. For a van- 
ishing potential, this equation has only a decaying and 
a constant solution; thus, the value of the field is fixed 
in time in each causally disconnected region of the early 
Universe. This precludes the scalar-field perturbations 
from having any correlation with perturbations in the 
matter / radiation density. This is manifest in the absence 
of any source term in the perturbed equation of motion 
for the scalar field (compare to the Fourier transform of 
the full equation, Eq. (jA2p . after taking (ci^/dr) — 0, 
and V = 0), 



(4) 



A solution to Eq. Q is a transfer function Tfc(r) oc 
ji{kT)/ (kr), which describes the conformal-time evolu- 
tion of a given Fourier mode of wavenumber k during 
matter domination. 

The angular power spectrum C£" for the rotation an- 
gle is then 



M 



mi 



2tt'^ 



^ P^k) [jL{kAT)Tkims)]' 



M 



(5) 



Here Ar is the conformal-time difference between last 
scattering and today, and tiss is the conformal time at 



^ Wc imagine that some mechanism has nuUified the quantum- 
gravity effects that generically break global symmetries l26l . 

^ It is also imaginable that a white-noise spectrum of (j> fluctuations 
is imprinted by some post-inflation phase transition, but we will 
not consider that scenario here. 
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the LSS. For large angular scales, L < 100, the transfer 
function evaluates to T]^{Tigg) ~ 1, in which case 



C 



27rL(i + 1)' 



for L < 100. 



(6) 



The left-hand panel of Fig. [T] shows the result of a nu- 
merical calculation of C£" for this scale-invariant power 
spectrum. 

The mean-square rotation amplitude measured by a 
probe with angular resolution of ^ 1° is 



L=2 



332 {^) deg^. 



(7) 



Here, Wl{6) = exp [— L^0^/(161n2)] is a Gaussian win- 
dow function of full- width half- maximum (in radians). 
The best current constraint to the variance of the rota- 
tion angle, ((Aa)^) < (4°)^, from AGN data [13, places 
a bound jiHj/M < 0.2 to the combination of parame- 
ters that control the rotation-angle amplitude in this sce- 
nario. We may therefore write the rotation-angle power 
spectrum as 



0.015 al I ^[jL{kAT)Tk{Tis)f 



7.7 X 10-3al 
L{L + l) ■■ 



forL < 100. 



(8) 



1 /2 

Here, ^4 = {{Aaf) ' /A" parametrizes the amplitude of 
the rotation-angle power spectrum in units of the maxi- 
mum value currently allowed. 

The CMB temperature power spectrum is given by 



C 



TT 



k^dk [AT,L(fc)] P*(fc), 



(9) 



where P^{k) is the primordial power spectrum for the 
gravitational potential, and Ax^hik) is the transfer func- 
tion that quantifies the contribution of a density mode 
of wavenumber k to Cj^, and may be obtained from 
numerical Boltzmann codes [27j . 

As discussed above, scalar-field fluctuations are not 
sourced by the gravitational potentials for this V ~ 
model; similarly, energy-density fluctuations in the scalar 
fleld have only second-order corrections due to S(f), and 
so their effect on gravitational potentials is also small. In 
this case, the aT cross-correlation power spectrum van- 
ishes, Cl^ = 0. 



B. Quintessence 

In the second scenario, we suppose that is a 
quintessence fleld with a nonzero potential and homo- 
geneous component that undergoes time evolution. In 



this case, gravitational-potential perturbations directly 
source (and are also sourced by) scalar-fleld fluctuations, 
see Eq. (jA2| . A cross-correlation between CB-angle and 
CMB-temperature fluctuations is therefore inevitable, al- 
though its amplitude and detailed features depend on the 
speciflc potential V. 

Since every CMB photon that comes from a given di- 
rection n last scattered at the spacctime point in the 
direction n, when the Universe had some fixed temper- 
ature, the CB-rotation angle a(n) is determined by the 
value of (f> at that point of spacetime. In other words, the 
CB-angle anisotropies are determined by the scalar-field 
perturbations on surfaces of constant CMB temperature, 
or equivalently, on surfaces of constant synchronous- 
gauge time. 

We suppose that the initial value of is set by some 
post-inflationary physics so that the primordial perturba- 
tion to (j) is adiabatic. In this case, the synchronous-gauge 
scalar-fleld perturbation ((5^)syn is initially zero. How- 
ever, the scalar-field perturbation is sourced by the grav- 
itational potentials, as described by Eqs. (jAip or (|A2p . 
The synchronous-gauge scalar-field perturbation at the 
LSS is then approximately (sec Appendix [A]) . 

(')'/>)sy„,lss = - g ( ^^-^ I A/pi*, (10) 

where the equation-of-state parameter and the 
energy-density parameter are evaluated at recombi- 
nation. The primordial power spectrum for the gravita- 
tional potential, for large scales (small k) is given by 



25~k^ 



(11) 



where we have taken a scalar spectral index to be = 1 
for simplicity, and the curvature-perturbation amplitude 



is A^{ko) = 2.43(±0.11) X 10" 



To evolve the power 



spectrum from primordial to the LSS, we need to multiply 
it by transfer functions, which are a suppression factor 
for small scales (large fc's). The angular power spectrum 
for the CB-rotation angle in the quintessence model is 
then 



27 



$1^(1 + w^) 



//3M] 



PI 



fc2 dk 
27r2 



V M 



P^{k)[jL{kAT)Tk{Tiss)f. (12) 



For large scales, L < 100, we can approximate Tk{Tiss) — 
1, in which case we can again write the CB-rotation-angle 
power spectrum as in Eq. ([8]), but now with 



a4 ~ 6.7 X 10" 



' ^n^{l + w^){l3Mpi/M). (13) 



In other words, the aa power spectrum for the 
quintessence scenario will be similar to that for the 
massless-scalar-field scenario in the small-L limit where 
Tk{Tiss) can be approximated as a constant. 
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Figure 1: Shown are the power spectra for the cosmological- 
birefringence rotation angle C£" and its cross-correlation with 
the CMB temperature (logarithm of the absolute value), 
for a generic quintessence model in which the CB-angle fluc- 
tuations are due to scalar-field fluctuations at the LSS. The 
black solid curves are the theoretical prediction for (from top 
to bottom) a4 = 1, O.f, and 0.01, where a4 is the fluctuation 
amplitude for the CB angle in units of the maximum currently 
allowed amplitude [T^]. We also show the noise power spec- 
tra anticipated for SPIDER (red, dot-dashed), Planck (green, 
dashed), and CMBPol (blue, dotted). 



show how the cross-correlation with the temperature can 
be reconstructed. At the end, we evaluate the detectabil- 
ity of the CB-rotation with both the auto-correlation and 
the cross-correlation for the CB scenarios discussed in 
Section im 



A. Measuring the rotation angle 

Refs. [13, [2l| show how the CB-rotation-angle 
spherical-harmonic coefficients Ulm can be reconstructed 
from a full-sky CMB temperature/polarization map. 
While these coefficients can be obtained from EE, TE, 
TB, and EB cross-correlations, the best sensitivity will 
ultimately come from the EB cross-correlation. We there- 
fore restrict our attention to reconstruction of a(fi) from 
the EB power spectra. 

To begin, the E/B spherical-harmonic coefficients, 
^"m^ ^^'^ ^tm^ ' ^^'^ constructed from the full-sky map of 
the Stokes parameters, Q(n) and U{n), in the usual way 
[1^, d^l. Following Refs. [lo, [HI, the minimum- variance 
estimator for the rotation-angle spherical-harmonic coef- 
ficient is 



aLM = 



aa, noise \ ^ ^LM rr j^map nmap 



rnm' ^l' >l 



where 



F, 



L.BE 



BB, mcip y^EE , map ' 



(15) 



(16) 



However, in the quintessence model, there will also be 
a cross-correlation with the CMB temperature, since the 
CMB temperature is determined largely by the potential 
* at the LSS. From Eqs. dH) and dnj, we get 



47r n^{l + w^) ^Mpi 
' 3 V 6tt M 



dk 
27r2 



P^{k)AT,Lik)jL{kAT)Tk{Tis 



(14) 



The absolute value of this cross-correlation is also shown 
in Fig. [T] The passage through zero at i ^ 50 arises 
because of the relative contributions of the monopole and 
dipole contributions of the photon distribution function 
to At,l(/c). 



III. PROSPECTS FOR DETECTION 

In this Section, we first review the procedure presented 
in Refs. [20l - |22| for constructing the minimum- variance 
estimator for the CB-rotation angle from a CMB tem- 
perature/polarization map. We then extend this work to 



pL.BE _ 2Cjf^ 



and 



^Inil'ni' 



1 L I' 

2 0-2 



= F, 



L,BE 
VI ' 

(17) 



^ (-1)' 



(2Z + 1)(2Z' + 1)(2L+1) 



47r 



I L V 

-m A I m' 



(18) 



Here, the objects in parentheses are Wigner-3j symbols, 
and Wi is the window function defined in Section [Hi 
and O; are, respectively, power spectra 

for the E and B modes from the map (including instru- 
mental noise); i.e., 



^X.Y',map_^XX'|w|2 



a 



XX', 



(19) 



where XX' e {TT, EE, BB, ET, EB,TB}. The noise 
power spectra are 



C, 



.TT. noise 



47r/o,y(NET)^ 



^obs 



y^EE, noise /~iBB. noise n/^TT. noise 

(20) 



EB, 



c. 



TB, 
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where ^obs is the total observation time, f^^^ is the frac- 
tion of the sky surveyed (taken to be different from 1 only 
for SPIDER, where /^^y = 0.5), and NET is the noise- 
equivalent temperature. We assume no cross-correlation 
between the noises in polarization and temperature and 
apply the null assumption (no B modes in the signal), so 
there are no TB and EB correlations. The power spec- 
trum C;^^' thus contains only the contribution from 
instrumental noise. 

Under the null hypothesis of no rotation, the expec- 
tation value of the estimator in Eq. (|15p is zero, and 
its variance is the aa noise power spectrum as given in 
Ref. 



C" 



■aa, noise 



^ (2H- 
L W 



mi' 



l)(F^/^^)2 



^^^B B , ma,p^EE^ map 



.(21) 



{Im} = {LM}. There is therefore no correlation (un- 
der the null hypothesis) of Slm and T^m; i-e., there is 
no noise contribution to C£"^. Again, if a^M is approxi- 
mately Gaussian, then the variance with which can 
be measured is approximately that obtained assuming 
ctLM is Gaussian. To check the validity of this assump- 
tion for the purpose of calculating the sensitivity of future 
CMB experiments to the CB signal (see Section IIII Cp , 
we evaluate the full expression for this variance (without 
assuming Gaussianity of ctLM, see Eq. (|Bip ) and confirm 
that the numerical results agree up to a level of a few 
percent. Thus, for simplicity, and without any loss in ac- 
curacy, we can invoke analogy with the variance of C^^ 
(see, e.g., Ref. [H]) to get 



1 



/sky(2L+l) 



c 



aa, noise 



(25) 



If the polarization pattern at the LSS is a realization of 
a statistically isotropic field, then there are 2L + 1 statis- 
tically independent M modes for each L in ctLM- In this 
case, each M mode provides an independent estimator 
of the rotation power spectrum, C£". The minimum- 
variance estimator is then 



2L 



1 ^ 

— y 

+ 1 . 



\aLM\ 



(22) 



M=-L 



Each CtLM is a sum of products of Gaussian random vari- 
ables and is thus not a Gaussian random variable. How- 
ever, if the number of terms in the sum is large, the 
central-limit theorem holds, and Slai can be approxi- 
mated as Gaussian. In this case, the expression for the 
variance of C£" takes on the usual form. 



(Acr)' 



Uy{2L + 1) 



C" 



Q a, noise 



(23) 



where /sky is the sky-cut used in the analysis, taken to 
be 0.8 for Planck and CMBPol and 0.5 for SPIDER. 



B. Measurement of the rotation-temperature 
cross-correlation 



In analogy with the derivation in Ref. [29| of the esti- 
mator for C^^ , the estimator for is 



L 



2L+1 



J2 SLMiTZ7rwr:\ 

M=-L 



(24) 



where T^^m temperature spherical-harmonic coef- 

ficient obtained from the map. Under the null hypoth- 
esis, Tlm has no correlation with any B^ms, and it is 
correlated with Elj^i with the same L and M, but un- 
correlated with any other E^m- The estimator a^M de- 
pends on a large number of Ei^s but does not include 



C. Sensitivity to Detection: aT vs. aa 

We now return to our two models for CB which pre- 
dict that the rotation a is a realization of a random field 
with the power spectra C£" and presented in Fig.[T] 
Our aim here is to evaluate the smallest signal amplitude 
detectable by measurement of the rotation alone, as well 
as the smallest amplitude detectable by measurement of 
the rotation-temperature cross-correlation. 

We write the power spectra as Cff 



alC\ 



ta. fiducial 



and = a4C2"^'^'^^"^^ , where the fiducial model {a4 = 
1) is the quintessence model in Fig. [T] with the largest 
amplitude allowed by current rotation-angle constraints. 
The inverse- variance with which the amplitude a^ of the 
aa power spectrum can be obtained from the rotation- 
angle auto-correlation is [3l| 



V fdcr 



1 



L 



^aa, fiducial 



(26) 



Similarly, the inverse- variance with which the amplitude 
a4 of the aT power spectrum can be obtained from the 
cross-correlation of the rotation with the temperature is 



1 



[AaiY 




(27) 



From these relations, we can estimate the signal-to-noise 
ratio for measurement of a\ from the aa autocorrela- 
tion to be {S/N)aa = oi\/ \^A{a\)\ and a signal-to-noise 
for measurement of 0:4 from the aT cross-correlation to 



6 



be {S/N)aT = Q;4/(Aa4). We evaluate these expres- 
sions for our fiducial model (a^ = 1), for different in- 
strumental parameters in Section IIII Dl The smallest 
detectable at the 2a level from the cross-correlation and 
auto-correlation are then 2Aa4 and [2A(q;4)] respec- 
tively. 



Instrument 


e 


NET 


^obs 




(S/iV)cT 


SPIDER 


60 


3.1 


0.016 


9 


7 


Planck 


7.1 


62 


1.2 


11 


9 


CMBPol/EPIC 


5 


2.8 


4 


2 X lO'' 


1200 



Table I: Instrumental parameters from Refs. [2^ - [25t for 
the three experiments considered in this work: beamwidth 
9 (in arcminutes), noise-equivalent temperature (NET) (in 
/iK sec^'''^), and observation time tobs (in years). The last 
two columns list signal-to-noise ratios (S/N) for the CB-angle 
auto-correlation and its cross-correlation with the CMB tem- 
perature, for our fiducial quintessence model (04 = 1) shown 
in Fig. [1] Note that the signal-to-noise scales with the signal 
amplitude Q4 as {S/N)aa oc a| and {S/N)aT 04. 



D. Numerical Results 

We now present numerical results for the aa and aT 
noise power spectra and evaluate the largest possible 
signal-to-noise and the smallest detectable amplitude 
for three CMB polarization experiments: (i) SPIDER's 
150 GHz channel [H, (ii) Planck's 143 GHz channel [H, 
and (iii) CMBPol's (EPIC-2m) 150 GHz channel We 
obtain the CMB temperature-polarization power spectra 
from CMBFAST ^ using WMAP-7 cosmological pa- 
rameters (isj . The instrumental parameters we use arc 
listed in Table HI Fig. [T] shows the noise power spectra'^ 
^aa, noise ^aT.noisc^ C^^ , strictly Speaking, there 
is no instrumental-noise contribution, only the effective 
noise, arising from cosmic variance. Table |T] lists signal- 
to-noise ratios, assuming 04 = 1, for the auto- and cross- 
correlations, for these three experiments. We find that 
SPIDER and Planck may already have the sensitivity to 
detect not only the signal, but also its cross-correlation 
with with the temperature, in the best-case scenario of 
0:4 ~ 1, where the signal is just below the current detec- 
tion limit'''. In both cases, the sensitivity to the signal 



^ Note that there is a differe,nce in normahzation between the noise 
and the variance: C^^ = ^{2L + l)/2ACf^' , where 

XX' = {aa, aT}. It is customary to plot the noise power spec- 
tra, even though the variance enters the expressions for signal- 
to-noise. 

* Here we have assumed that the errors to the rotation-angle esti- 
mators are approximately Gaussian. However, if the signal is just 
detectable (for example, for experiments like SPIDER if 04 ~ 1), 
then this assumption may break down, and if so, the precise 
quantitative forecasts for the signal-to-noise may differ slightly. 



may be improved if the auto- and cross-correlations are 
measured in tandem. CMBPol should have sensitivity 
to a signal as small as 04 ~ 10~^, and a detection of 
the cross-correlation of very high signal-to-noise may be 
obtained with CMBPol if 04 ~ 1. 



IV. SUMMARY AND DISCUSSION 

If a quintessence field gives rise to cosmological bire- 
fringence, then a correlation between CB-rotation-angle 
fluctuations and CMB-temperature fluctuations is in- 
evitable. We calculated that cross-correlation assuming 
the initial quintessence perturbations are adiabatic. We 
also discussed, by way of contrast, a scenario in which 
the CB-inducing field is just a massless scalar field that 
has no correlation with primordial perturbations. 

We derived the minimum-variance estimator for the 
aT power spectrum that can be obtained from a CMB 
temperature-polarization map. We find that measure- 
ment of this cross-correlation may improve sensitivity to 
the CB signal in some cases where the signal would oth- 
erwise be only marginally detectable. We further show 
that a high signal-to-noise measurement of this cross- 
correlation is conceivable with forthcoming and future 
CMB experiments if the rotation-angle power-spectrum 
amplitude is near its current upper limit. Measurement 
of this cross-correlation may thus provide another empir- 
ical handle with which to discover new physics indicated 
by cosmological birefringence. 

We have restricted our attention to the EB estimator 
for the rotation angle, as it is expected to provide the 
best sensitivity. However, there may be some improve- 
ment, though probably small, with the inclusion of the 
TE, TB, and EE estimators for the rotation. We leave 
this calculation for future work. Likewise, we have left 
more careful investigation of the impacts of partial-sky 
analysis, foregrounds, uneven noise, as well as the de- 
tails of distinction of these signals from weak lensing, for 
future work. 

We have refrained from discussing details of the 
quintessence model here, as the angular dependence of 
the CB power spectra at superhorizon scales at the time 
of recombination, L < 100, is insensitive to these de- 
tails. The dependence of the amplitudes of the aa and 
aT power spectra is given in terms of the quintessence 
parameters fi^ and w at the LSS by Eqs. and ([14]). 
However, if the quintessence field couples to the pseu- 
doscalar of electromagnetism, it is natural to expect it to 
be a pseudo-Nambu-Goldstone field, and if so, then its 
potential should be V{(j)) oc [1 — cos(0//)]. In this case, 
the quintessence field (f) is frozen at early times leading 
to spatial variations in a that are unobservably small. 
In this case, though, additional fluctuations in a may be 
produced during the epoch of reionization [3^ . 

For the massless scalar field, the uniform CB-rotation 
angle is expected to be zero, and so a search for the fluctu- 
ations is essential to detect the signal. For quintessence, 
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however, the uniform rotation is expected to be nonzero 
and generically quite a bit larger than the fluctuations, 
which, given the current best constraint may imply a rel- 
atively small amplitude of the fluctuations power spec- 
trum. However, CMBPol may be sensitive to a fluctua- 
tion amplitude as small as ~ 10^^ of the current upper 
limit to the uniform rotation, which, if detected, would 
help distinguish between different CB scenarios. More- 
over, the fluctuation amplitude in the quintessence sce- 
nario could be larger than a measured uniform-rotation 
angle. This could occur if, for example, the uniform- 
rotation angle (which can only be recovered mod tt) hap- 
pens to be close to an integer multiple of tt. It will 
be interesting, with forthcoming precise CMB maps, to 
address these questions empirically rather than through 
theoretical speculation. 
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Appendix A: Quintessence Perturbations 

As discussed in the paper, the CB-angle fluctuation is 
determined by the synchronous-gauge scalar-field fiuctu- 
ation {S(j))syn at the LSS. To obtain this fiuctuation, we 
start from adiabatic initial conditions, and then evolve 
the scalar-field-perturbation equation of motion forward 
in time, from the early radiation-dominated epoch to the 
LSS. The equation of motion is 



2H6<p + a^V"d<p - V 



1. 



-h(p, 



(Al) 



in the synchronous gauge, and 

S^+2n6(j)+a^V"S(l)~\7'^S(j) = (/)(3$-|-*)-2aV*, (A2) 

in the conformal-Newtonian/longitudinal gauge. (See 
Ref. [l^l for definitions of the metric variables $, 77, 
and h.) Adiabatic initial conditions require that the per- 
turbations of the scalar field vanish at early times. How- 
ever, the subsequent evolution of the scalar field is not 
adiabatic, meaning that ((5(?!))syn.iss docs not necessarily 
vanish at the LSS, even though all the matter and radia- 
tion perturbations do. All of our numerical integrations 
that give the power spectra presented in the Figure are 
done in the s ync hronous gauge, using a modified version 
of CMBFast 123. 

For the numerical work, we assume a quintessence po- 
tential of the PNGB form V{(j)) = m'^{l - cos(^//), as 



expected if (f) is an axion-like field. We take an ini- 
tial value (j), m, and / so that fi^ = 0.7 today and 
the density-weighted average equation-of-state parame- 
ter (w) ~ -0.95, which gives A(j> = 0.045 Mpi for the 
change in the scalar field between decoupling and today. 
However, the numerical results presented in the Figure 
will be similar for any quintessence potential that has 
— > — 1 at early times. 
The numerical results can be largely reproduced 
with the analytic approximation for (5</>)syn,issi given in 
Eq. (jlOp . which we now derive. We now work in the 
conformal-Newtonian/longitudinal gauge and we make 
the approximation that decoupling takes place well into 
matter domination; we assume that most of the growth in 
perturbations happens during this epoch. For ^ — 1, 
the V" term in Eq. (|A2p is negligible. Additionally, in the 
super horizon limit, valid for multipoles L < 100, we can 
neglect the spatial-gradient term. The simplified equa- 
tion of motion is then. 



2n5(j) 



-2a^Y'^. 



(A3) 



Aside from the homogeneous solutions that are either 
constant or decaying, it also has an inhomogeneous solu- 
tion that grows as 



-aW'«'/27, 



(A4) 



during matter domination. The potential derivative V 
can be expressed, using the quintessence slow-roll approx- 
imation, from 



a'V' 



-3-H0. 



(A5) 



Also, 



!i? = a^p^{\^w^), (A6) 
with /90 = Jl^pc and pc ~ 3i/^Mp;/(87r). We then find 



1/2 



Mpi^ 



(A7) 



The result in Eq. (|10p is then obtained by going back to 
the synchronous gauge, using the gauge-transformation 
equations [33^ 



m 



con - Otq 
con 



(A8) 
(A9) 



after noting that a ~ (2/3)^'/H, during matter domina- 
tion. 

We derive the initial conditions in the conformal- 
Newtonian/longitudinal gauge, for the sake of complete- 
ness, which can then be used to evolve Eq. (jA2p . To ob- 
tain the initial conditions in this gauge, we use Eq. (jA9p , 
where a = (1/2)"^ /Ti. during radiation domination. At 
early times, deep in the radiation era, we can set the frac- 
tional energy- density perturbation in the radiation field 
to 5r = —2^ [3I], and assume that the equation-of-state 
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parameter ^ — 1, and changes slowly with time. Fur- 
thermore, the pressure and energy density of the scalar 
are given as 



1 
2a 



1 

2^ 



0' + V{cb), 



(AlO) 



1 



and the perturbation in the energy density is 

iSp4.)con = \4>U + V'{4>)5^ - ^02$, (All) 

while the adiabatic initial conditions require that the en- 
tropy density perturbation vanishes at early times, so 
that 



S 



5pr 



Pr +Pr 



(A12) 



Combining these assumptions into the gauge- 
transformation equations, along with the observation 
that 6(j) vanishes in the synchronous gauge, we get the 
initial conditions for the scalar-field perturbations in the 
conformal-Newtonian /longitudinal gauge , 



These initial conditions can also be derived by requiring 
that S and S vanish at early times. 

Appendix B: Full expression for the variance of aT 
cross-correlation 

If we do not assume a^M is a Gaussian, then the 
full expression for the variance of its cross-correlation 
with the CMB temperature becomes a 6-point correla- 
tion function. After applying Wick's theorem and taking 
into account the properties of the Wigner 3j symbols to 
simplify the terms, the full expression becomes 



aT\2 



{Acr ) 



aa^noisc \ 2 j^B B . noise 



4:TtWI 



{21 + 1)^ fT^L\2r'EE, map^TT, map 



II 



(2^ + i)(2r + i) 

(2L + 1) 



L\2f-<TT, map qEE, map 



(Bl) 



2n ' 



(A13) 
(A14) 



[1] S. M. Carroll, G. B. Field and R. Jackiw, Phys. Rev. D 

41, 1231 (1990). 
[2] B. Ratra and P. J. E. Peebles, Phys. Rev. D 37, 
3406 (1988); C. Wetterich, Astron. Astrophys. 301, 
321 (1995) (arXiv:hep-th /9408025 ; K. Coble, S. Dodel- 
son and J. A. Frieman, Phys. Rev. D 55, 1851 (1997) 
[arXiv:astro-ph/9608122 ; M. S. Turner and M. J. White, 
Phys. Rev. D 56, 4439 (1997) arXiv:astro-ph/9701138 ; 
R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. 
Lett. 80, 1582 (1998) arXiv:a stro-ph/9708069] . 
[3] R. R. Caldwell and M. K amionkowski , Ann. Rev. 
Nucl. Part. Sci. 59, 397 (2009) [arXiv:0903. 0866 [astro- 
ph.CO]]; E. J. Copeland, M. Sami and S. Tsujikawa, Int. 
J. Mod. Phys. D 15, 1753 (2006) arXiv:hep-th/0603057 ; 
A. Silvestri and M. Trodden, Rept. Prog. Phys. 72, 
096901 (2009) arXiv:0904.0024 [astro-ph.CO]]. 
[4] A. G. Riess et al. [Supernova Search Team 
Collaboration], Astron . J. 116, 1009 (1998) 
'arXiv:astro-ph/9805201 . 
[5] S. Perlmutter et al. [Supernova Cosmology Project 
Collaboration], Astrophys. J. 517, 565 (1999) 
arXiv:astr o-ph/9812133.. 



[6] A. Albrecht et al, arXiv:astro-ph/060959l| J. Frieman, 
M. S. Turner and D. Huterer, Ann. Rev. Astron. As- 
trophys. 46, 385 (2008) [arXiv:0803. 09821 [astro-ph]]; 



E . V . Under, Rept. Prog. Phys. 71, 056901 (2008) 
arXwTOSO 1.2968 [astro-ph]]. 

[7] S. M. Carroll, Phys. Rev. Lett. 81, 3067 (1998) 
[arXiv:astro-ph/9 806099 . 

[8] J. A. Frieman, C. T. Hill, A Stebbins and I. W aga, Phys. 
Rev. Lett. 75, 2077 (1995) [arXiv:astro-ph79 505060 . 

[9] M. Pospelov, A. R itz, C. Sk ordis, Phys. Rev. Lett. 
103, 051302 (2009) |arXiv:080 8.0673 [astro-ph]]; M. Li 
and X. Zhang, Phys. Rev. D 78, 103516 (2008) 
, arXiv:0810.0403 [astro-ph]] . 
[10] S. Gardner, Phys. Rev . Lett. 100, 041303 (2008) 

'arXiv:astro -ph/0611684| . 
[11] B. Nodland and J. P. Ralston, Phys . Rev. Lett. 78, 
3043 (1997) arXiv:astro-ph/9704196;; T. J. Loredo, 
E. E. Flanagan and I. M. Wasserman, Phys. Rev. D 
56, 7507 (1997) [arXiv:astro-ph/9706258| ; D. J. Eisen- 
stein an d E. F. Bunn, Phys. Rev. Lett. 79, 1957 
(1997) |arXiv:astro-ph/9704247| ; J. F. C. War- 
die, R. A. Perley and M. H. Cohen, Phys. Rev . 
Lett. 79, 1801 (1997) ar Xiv:astro-ph/9705142] ; 
J. P. Leahy, arXiv:astro-ph/9704285' S. M. Carroll 
and G. B. Field, Phys. Rev. Lett. 79, 2394 (1997) 



arXiv:astro-ph /9704263] ; A. Cimatti, S. di Serego 



Alighieri, G. B. Field and R. A. E. Fosbury, Astro- 
phys. J. 422, 562 (1994); S. d. S. Alighieri, F. FineUi 



9 



[12] 
[13] 

[14] 



and M. Galaverni, Astrophys. J. 715, 33 (2010) 
[arXiv: 1003.4823 [astro-ph.CO]]; P. P. Kronberg, 
C. C. Dyer, and H.-J. Roser, Astrophys. J. 472, 115 
(1996); F. Fin elli and M . Galaverni, Phys. Rev. D 79, 
063002 (2009) l' arXiv:080 2.4210 [astro-ph]]. 
M. Kamionkowski, Phys. Rev. D 82, 047302 (2010) 
arXiv: 1004.3544 [astro-ph.CO]]. 

A. Lue, L. M. Wang and M. Kamionkowski, Phys. 
Rev. Lett. 8 3, 1506 (1999) ^arXiv:astro-ph/9812088, ; 
N. F. Lepora, [arXiv:gr-qc/9812077l 



[15] 
[16] 

[17] 

[18] 
[19] 
[20] 
[21] 



B. Feng, M. Li, J. Q. Xia, X. Chen and X. Zhang, Phys. 
Rev. Lett. 96, 221302 (2006) arXiv:astro-ph/0601095 ; 
P. Cabella, P. Natoh and J. Silk, Phys. Rev. D 76, 
123014 (2007) arXiv:Q7Q5.0 810 [astro-ph]]; T. Kahni- 
ashvili, R. Durrer and Y. Maravin, Phys. Rev. D 78, 
123009 (2008) |arXiv:0807.2"593 [astro-ph]]; J. Q. Xia, 
H. Li, G. B. Zhao and X. Zhang, Astrophys. J. 679, 
L61 (2008) >rXiv:0803.2350 [astro-ph]] ; J. Q. Xia, H. Li, 
X. 1. Wang and X. m. Zhang, Astron. Astrophys. 483, 
715 (2008) arXiv:0710.3325 [hep-ph]]; L. Pagano et ai, 
Phys. Rev. D 80, 043522 (2009) arXiv:0905.1651 [astro- 
ph.CO]]. 

E. Komatsu et al, arXiv:1001.4538' [astro-ph.CO]. 

H. C. Chiang et al, Astrophys. J. 711, 1123 (2010) 

iarXiv:0906.1181 [astro-ph.CO]]. 

J. Q. Xia, H. Li and X. Zhang, a rXiv:0908.1876l [astro- 
ph.CO]. 

E. Y. S. Wu et al. [QU aD Collaboration] , Phys. Rev. 
Lett. 102, 161302 (2009) |arXiv:0811.06T8l [astro-ph]]. 
M. L. Brown et al. [QUaD collaboration], Astrophys. J. 
705, 978 (2009) arXiv:0906.1003 [astro-ph.CO]]. 
M. Kamionkowski, Phys. Rev. Lett. 102, 111302 (2009) 
arXiv:0810.1286 [astro-ph]]. 

V. Gluscevic, M. Kamionkowski and A. Cooray, Phys. 
Rev. D 80, 023510 (2009) [arXiv:0905.1687l [astro- 



[22] 



[23] 
[24] 
[25] 

[26] 



[27] 

[28] 
[29] 



[30] 
[31] 

[32] 
[33] 



ph.CO]]. 

A. P. S. Yadav, R. Biswas, M. Su and M. Zaldarriaga, 
Phys. Rev. D 79, 123009 (2009) |arXiv:0902.4466l [astro- 
ph.CO]]. 

B. P. Crih et aL,larXiv:0807. 1548 1 [astro-ph] . 
'http : //www . rssd . esa . int /PLANCK 

J. Bock et al, ^rXiv:0805.4207l [astro-ph]; 



http://cmbpol.uchicago.edu/ 



M. Kamionkowski and J. March-Russell, Phys. Lett. 
B 282, 137 (1992) [arXivihep-th /9202003| ; R. Holman, 
S. D. H. Hsu, T. W. Kephart, E. W. Kolb, R. Watkins 
and L. M. Widrow, Phys. Lett. B 282, 132 (1992) 
arXiv:hep-ph/9203206 . 

U. Seljak and M. Zaldarriaga, Astrophys. J. 469, 437 



(1996) [arXiv:astrcKph/9603033| 
D. Larson et aL. larXiv:1001.4635l [astro-ph.CO]. 
M. Kamionkowski, A. Kosowsky and A. Stebbins, 
Phys. Rev. D 55, 7368 (1997) 'arXiv:astro-ph/9611125|; 
M. Kamionkowski, A. Kosowsky and A. Stebbins, Phys. 
Rev. Lett. 78, 2058 (1997) |arXIvrastro^ph/9609132 . 
M. Zaldarriaga and U. Seljak, Phys. Rev. D 55, 
1830 (1997) |arXiv:astrcHph/9609170| ; U. Seljak and 
M. Zaldarriaga, Phys. Rev. Lett. 78, 2054 (1997) 
arXiv;astro-ph/9609169; . 

G. Jungman, M. Kamionkowski, A. Kosowsky 
and D. N. Spergel, Phys. Rev. D 54, 1332 
(1996) 



arXiv:astro-ph/9512139] ; 



G. Jungman, 

M. Kamionkowski, A. Kosowsky and D. N. Spergel, Phys. 
Rev. Lett. 76, 1007 (1996) [arXiv:astro^ h/9507080 . 
R. R. Caldwell, V. Gluscevic, and M. Kamionkowski, in 
preparation. 

C. P. Ma and E. Bertschinger, Astrophys. J. 455, 7 
(1995) jarXiv:astro-ph/9506072, . 



